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Abstract 

> 

\ The difFractive open charm production is computed in perturbative QCD formahsm 

■ and in the Regge approach. The results are compared with recent data on charm 

diffr active structure function measured at DESY-HERA. Our results demonstrate 
' that this observable can be useful to discriminate the QCD dynamics. 

o 



1 Introduction 



The study of electroproduction at small x has lead to the improvement of 
our understanding of QCD dynamics at the interface of perturbative and non- 
perturbative physics. In particular, the discovery of diffractive events in this 
process at HERA has triggered a large amount of experimental and theoretical 
work and greatly increased our knowledge of the physics of diffraction (For 
recent reviews see Refs. [1,2,3]). Diffractive processes in deep-inelastic scat- 
tering (DIS) are of particular interest, because the hard photon in the initial 
state gives rise to the hope that, at least in part, the scattering amplitude can 
be calculated in pQCD. Moreover, DIS exhibits the nice feature of having a 
colorless particle, the virtual photon, in the initial state. The main theoretical 
interest in diffraction is centered around the interplay between the soft and 
hard physics. Hard physics is associated with the well established parton pic- 
ture and perturbative QCD, and is applicable to processes for which a large 
scale is present. Soft dynamics on the other hand, linked for example with 
the total cross section of hadron scattering, is described by nonperturbative 
aspects of QCD. The ability to separate clearly the regimes dominated by soft 
and hard processes is essential in exploring QCD at both quantitative and 
qualitative level. 
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Recently, we have proposed the analyzes of the slope of diffractive structure 
function as a potential observable to disentangle the leading dynamics at ep 
diffractive processes [4,5]. The predictions for the behavior of this quantity are 
strongly dependent of the QCD dynamics dominant in the kinematical region 
(For a recent discussion see Ref. [6]). Similarly, the study of the diffractive 
final state can lead to further progress in the direction of obtaining a coherent 
picture of the diffraction. In particular, charm production looks promising in 
this respect, as predictions for this process widely differ among several models 
[7,8,9,10,11]. Recently, the ZEUS collaboration has presented its results for 
the measurement of the open-charm contribution to the diffractive proton 
structure function [12]. Consequently, a more detailed analyzes of the models 
and a comparison between their predictions and the experimental data is on 
time. Here the diffractive open charm production is computed in perturbative 
QCD formalism and in the Regge approach. As these models are based on 
very distinct assumptions, it allows shed light into the leading dynamics at ep 
diffractive processes. 

This paper is organized as follows. In the next section, we summarize the 
main formulas for computation of the open charm diffractive structure func- 
tion. One presents it in the transverse momentum representation in the per- 
turbative QCD approach. Moreover, the diffractive production of open charm 
is calculated in a Regge inspired approach, where charm is produced by bo- 
son gluon fusion and which directly depends on the gluon distribution of the 
Pomeron. In the last section, we compare both approaches with current exper- 
imental measurements from DESY-HERA collider and present our discussions 
and conclusion. 



2 Diffractive production of open chcirni in deep inelastic scattering 



Before introducing the main expressions needed to our calculation, let's intro- 
duce the kinematical definitions in diffractive DIS (DDIS) 7* (g) + p (P) 
X (Mx) p (-P') with X being the diffractive final state. The kinematics is de- 
fined as follows. 



^ 2P.g' q.P ' ^~ 2g.(P-P0 ^ + 

where g, P and P' are the four-momenta of the virtual boson, the incident 

proton and the remnant colorless final state, respectively. The invariant mass 
of the diffractive final state is labeled Mx- The variable x is the momentum 
fraction of the proton carried by the partons (quarks or gluons), the Bjorken 
variable, and by definition x — (5xip. As usual, — —(f is the photon 
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virtuality. 



At high energies, xjp may be interpreted as the fraction of the proton four- 
momentum carried by the diffr active exchange, the colorless Pomeron. The /3 
variable is the fraction of the four-momentum of the diffractive exchange car- 
ried by the parton interacting with the virtual boson. The diffractive structure 
function is defined in analogy with the decomposition of the unpolarized total 
ep cross section as. 



Since the first observation of diffractive DIS at HERA, several attempts have 
been made to compare the data with the Regge and QCD-based models 
[14,15,16] (See also [17,18,19]). In general, these models provide a reasonable 
description of the present data on the diffractive structure function F^^^\ 
although based on quite distinct frameworks. Furthermore, the QCD factor- 
ization theorem has been proven to be valid for F^^^^^ [20], with the imme- 
diate consequence that the DGLAP evolution equations [21] should describe 
the scaling violations observed in this observable. However, there are no con- 
straints on the gluon momentum distribution since the momentum sum rule 
does not formally apply. 

Here, we study in detail the predictions for the charm component of the diffrac- 
tive structure function, jp^i^)^^^^^ which is directly sensitive to the gluonic 
content of the pomeron, considering two distinct approaches: i) a Regge in- 
spired model [13,14], where the diffractive production is dominated by a non- 
perturbativc Pomeron, and the diffractive structure function is obtained us- 
ing the Ingelman-Schlein ansatz [22]. ii) a pQCD approach [15,16] where the 
diffractive process is modeled as the scattering of the photon Fock states with 
the proton through a gluon ladder exchange (in the proton rest frame) . Below 
we present a brief review of the main assumptions of these models. 

In the perturbative QCD framework, there are successful analysis describing 
the diffractive structure function [15,16]. The underlying physical picture is 
that, in the proton rest frame, diffractive DIS is described by the interaction 
of the photon Fock states {qq and qqg configurations) with the proton through 
a Pomeron exchange, modeled as a two hard gluon exchange. The correspond- 
ing structure function contains the contribution of qq production to both the 
longitudinal and the transverse polarization of the incoming photon and of 
the production of qqg final states from transverse photons. The basic elements 
of this approach are the photon light-cone wave function and the noninte- 
grated gluon distribution (or dipole cross section in the dipole formalism). 
For elementary quark-antiquark final state, the wave functions depend on the 
helicities of the photon and of the (anti) quark. For the qqg system one con- 
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siders a gluon dipolc, where the pair forms an effective gluon state associated 
in color to the emitted gluon and only the transverse photon polarization is 
important. The interaction with the proton target is modeled by two gluon 
exchange, where they couple in all possible combinations to the dipole. Then 
the diffractive structure function can be written as [15,16] 



(3) 



where D\E' is a combination of the concerned wave functions, It is the trans- 
verse momentum of the exchanged gluons. The function T{xjp. if) defines the 
Pomeron amplitude (nonintegrated gluon distribution) and contains all the de- 
tails concerning the coupling of the i-channel gluons to the proton. Integrating 
it over one obtains the conventional coUinear gluon distribution. 

Concerning diffractive open charm production, the exclusive cc-pair arises 
from the dissociation of longitudinally and transversely polarized photons, 
as well as the production of the ccg^-state. The diffractive structure functions 
for 7*p — > ccp are given by [7,8,1], 



xipF^Jxjp,P,Q') 



D (1 - py J 



A8B, 



Af + ml 



dk^ 



r dk 
3SbQ2 J 



+ 



k^ 



Pyji- Apkyq- 



>,(4) 

r(5) 



where the upper limit in the integration on the quark loop is constrained by 
the 0-function. The parameter Bd is the diffractive slope, which arises by 
assuming a simple exponential form for the \t\ dependence to the process (one 
uses Ed — G GeV~^ in the following). The integrals It,l on gluon transverse 
momentum are defined as. 



f d£'^ 



-(1-2/3) fc2 + 2m; 



fc2 



^(£2 + p)2_4£2^2j 



where the two-body kinematical relation has a mass term (in relation to the 



light quarks dipoles) and reads as A;^ = • The allowed range on /3 is 
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different from the light dipole case as the difFractive mass Mx has a lower 
limit defined by M"^ > 4m^. For the energy scale entering in the strong 
coupling we will use the prescription /^^ = 4m^. 

In order to compute the contribution of the ccg component, one makes use 
of the difFractive factorization property [1]. The difFractive gluon distribution 
g^{P) will be convoluted with the corresponding charm-coefficient function 

where the lower limit in the z integration is weighted by a = 1 + Arn^/Q^ and 
the coefficient function is given by, 



Cg (C, r) = [(' + (1 - 0' + 4C(1 - 3C)r - SCV^] In 

+ £[-l + 8C(l-C)-4C(l-C)r] , (7) 

with e, the centre-of-mass velocity of the charm quark or antiquark, given by 
£ = ^l-(4rC/l-C)- 

For the difFractive gluon distribution, we use the momentum representation, 
which reads as [1] 



X 



^2 ^ (1 _ ^ :f _ [(1 - - ft? + 2(3(1 - (3)k' 



+ A;2)2 - 4(1 - /?) 



2 

(8) 



The upper limit oF the /cj-integration is fixed by the condition Mx = Q"^ ^^ > 

{kt + sJkl + Amlf, which implies kj < 

The difFractive gluon distribution obtained above depends directly on the un- 
integrated gluon Function. Concerning the behavior on /5, an expansion in 
powers oF fjk^ [1] produces ^^{1- (3f{l + 2(3f J dk^ [xipg{xip, k'^)f/ki, 
where g{xip, Q'^) is the coUinear gluon distribution. ThereFore, the diFFractive 
gluon distribution has a singular behavior at /3 — > and vanishes at — > 1. 

In order to perForm Further numerical analysis, we will use the unintegrated 
gluon Function giving by the saturation model, which has a simple analytical 
Form [16]. It reads as. 
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Fig. 1. The open charm diffractive structure function and its three components 
plotted versus /3 for a fixed xp = 0.004 and = 4 GeV^. 

where QLtl^) — f~) ' '^sat is the saturation scale and one has used the 
parameters for the 4-flavor fit. Accordingly, for the computation of the qqg 
contribution, Eq. (9) has been properly rescaled concerning the color charge 
[23]. 

In Fig. 1 we present how the three contributions, ccg, cc from transversely and 
longitudinally polarized photons, contribute for the /3-spectrum of ^^^^"^ _ 
At small-/? we have that the ccg component dominates, which implies that the 
fraction of charm in this regime is predicted to be the same as expected in 
inclusive charm production 25%). The above result agree with the theoret- 
ical expectations [1] . Since the mass of the quark sets a limit on the size of the 
cc dipole, it becomes color transparent and one expects a strong suppression 
for this configuration. On the other hand, the effective gluon dipole, associated 
with ccg production is not restricted in size. 

Concerning the Regge inspired approaches, diffraction dissociation of virtual 
photons furnishes the details on the nature of the Pomeron and on its partonic 
structure. As a first investigation, we follows the Capella-Kaidalov-Merino- 
Tran Thanh Van (CKMT) model to diffractive DIS based on Regge theory 
[13,14] and the Ingelman-Schlein ansatz, which is based on the intuitive pic- 
ture of a Pomeron flux associated with the proton beam and on the conven- 
tional partonic description of the Pomeron-photon collision. In this case, deep 
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inelastic diffractive scattering proceeds in two steps (the Regge factorization): 
first a Pomeron is emitted from the proton and then the virtual photon is 
absorbed by a constituent of the Pomeron, in the same way as the partonic 
structure of the hadrons. In the CKMT model the structure function of the 
Pomeron, Fip{P,Q^), is associated to the deuteron structure function. The 
Pomeron is considered as a Regge pole with a trajectory a]p{t) determined 
from soft processes, in which absorptive corrections (Regge cuts) are taken 
into account. The diffractive contribution to DIS is written in the factorized 
form, 

F,''^'\xjp,P,Q\t) = f{xjp,t)Fjp{P, Q') , (10) 

where the first factor represents the pomeron fiux from the proton and can be 
written as 



fM^^^x^-^'^^^^ (11) 



where g^p{t) = 5^(0) exp(C t) is the Pomeron-proton coupling, with [(/^(O)]^ = 
23 mb and C = 2.2 GeV~^ [13,14]. The Regge factorization implies that the 
X]p dependence is completely separated from the (3 dependence, with the be- 
havior in xip determined only by the fiux factor. The value of a]p{t) in the fiux 
is given by 



aip{t) ^ 1 + A{Qlff) + a' t , (12) 
where a' = 0.25 GeV^^ and 



^«^'-^(»>0+^) ■ '''' 

with A(0) = 0.09663, do = 1.9533 and di = 1.1606 [24]. The dependence of 
the effective Pomeron intercept is one of the main feature of the CKMT model. 
It was argued in the Refs. [13,14] that this is due to the fact that the size of 
the absorptive corrections decreases when increases. This parameterization 
gives a good description of all existing data on 7*p total cross section in the 
region < 10 GeV^ [24]. At larger Q"^, effects due to QCD evolution become 
important. The scale Qlff is a priori not known. From a theoretical point of 
view, values for A(Qe//) between 0.13 and 0.24 are possible, corresponding 
to the effective Pomeron intercept without eikonal-type corrections and the 
"bare" value, respectively. Both values are not excluded by the recent fit for 
the data which assumes in addition to the Pomeron exchange, the contribution 
of a subleading reggeon trajectory [18,19]. Integrating Eq. (10) over t, F^^^^ 
can be put in the factorized form 
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F.'^^'^xjp, P, Q') = fixjp) Fjp iP, Q') , 
where f{x]p) is the t-integrated pomeron flux 



(14) 



oo 

f{xip)^ j d\t\f{xjp,t) . (15) 



It must be stressed that since the Pomeron is not a particle the separation of 
the flux factor from the photon-Pomeron cross section is quite arbitrary, and 
therefore the normalization of the flux is ambiguous. 

The second factor in Eq. (10) is the pomeron structure function Fp and is 
proportional to the virtual photon-pomeron cross cross section. In the CKMT 

approach, Fip{(3, Q^) is determined using Regge factorization and the values 
of the triple Regge couplings determined from soft diffraction data. Namely, 
the Pomeron structure function is obtained from F2, or more precisely from 
the combination F2 — |(F2 +F^), by replacing the Reggeon-proton couplings 
by the corresponding triple reggeon couphngs (see Ref. [13] for details). The 
following paramctrization of the deuteron structure function F2 at moderate 
values of (and small-x) , based on Regge theory, was introduced. 



Fi{x,Q')=Ax-^(<^'\l-xr^<^'^^' (^^j , (16) 

where 1 + A{Q'^) is the Pomeron intercept, which depends on the photon vir- 
tuality. The Pomeron structure function is identical to F2 except for a simple 
changes in its parameters: F^.(/?, Q^) ^ F2 {x ^ (5\ A ^ eA, n{Q^) n{Q'^) - 2). 
The value of e in Fjp is obtained from conventional triple reggeon fits to high 
mass single diffraction dissociation for soft hadronic processes. The remaining 
parameters are given in Refs. [13,14]. In the CKMT approach, the gluon dis- 
tribution of the Pomeron can be obtained for low /3 in a similar way as for the 
quarks discussed above. It is written as. 



P9IP {f3, Q') ^^C,P~^^^'^ (1 - , (17) 

where Ug is a free parameter and e'^ = jp{t) / (fdd ~ 0.07, with r^jp and 
(/^^ being the couplings of the Pomeron to the Pomeron and to the deuteron, 
respectively. The distribution is singular towards /3 — > due to the powerlike 
behavior driven by the Pomeron exchange. In Ref. [11], where charm diffractive 
production was computed, Ug takes values between and -1 in order to produce 
a normalizablc distribution, which implies for n^, < a singular behavior also 
at /5 = 1. As a good description of the dependence of the HERA data is 
achieved with Ug — Q and in view that a singular behavior for large /? is not 
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observed in Regge-like fitting procedures to DDIS data, we assume this value 
in our analyzes. 

In the Regge based approaches, the massive charm contribution arises from 
photon-gluon fusion. The diffractive structure function is ^2^*'^'' = f{xip) x 
Fp{(3, Q^), where f{xjp) is given by Eq. (15) and the charmed contribution to 
the diffractive Pomeron structure function, Fp{(3, Q'^), is given by folding the 
gluon distribution Eq. (17) in Eq. (6) [11]. The factorization scale is assumed 
equal to 4m^. In the general case, the scahng violations of Pomeron struc- 
ture function should be considered. However, as the Q^-range of the HERA 
data considered here is either small and the parameters in (17) have been ob- 
tained for = 5 GeV^, in a first approximation we disregard the logarithmic 
dependence on Q^, which is given by QCD-evolution. 

Another possible approach is the QCD analysis of the diffractive structure 
function in terms of both Regge phenomenology and perturbative QCD evo- 
lution as made in Ref. [18]. In this case the parton distributions of the Pomeron 
are derived from QCD fits of diffrative deep inelastic scattering cross sections 
determined at HERA. In particular, the diffractive structure functions is given 

by 



F,''^'\q',(3,xjp) = fjp/,{xjp)F^{Q',(3) + fmUxip)F^{Q\P) , (18) 

where F^ can be interpreted as the Pomeron structure function and F^ as 
an effective Reggeon structure function, with the restriction that it takes into 
account various secondary Regge contributions which can hardly be separated. 
The Pomeron and Reggeon fiuxes are assumed to follow a Regge behavior with 
linear trajectories aip^{t) — a]p^{0) + a'jp^t, such that 



tcut 

where \tmin\ is the minimum kinematically allowed value oi \t\ and tcut — ~^ 
GeV^ is the limit of the measurement. The gluon distribution for the Pomeron 
is parameterized in terms of non-perturbative input distributions at Qo = 3 
GeV^ as foUows 



f3G{(3,Q' = Ql) 



(20) 



and similarly for the quark fiavor singlet distribution. The Pj{f]) is the j'*'* 
member in a set of Chebyshev polynomials, which are chosen such that Pi — 1, 
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Fig. 2. The open charm diffractive structure function F2 as a function of (5 

(data from ZEUS Collaboration [12]). The solid lines correspond to the perturbative 
QCD calculation, whereas the other lines represents the results from the Regge 
approach. The long-dashed curves stand for the CKMT Pomeron and the dashed 
ones for a fixed Pomeron intercept (see text). The dot-dashed curves represents the 
result using a Regge/QCD analysis from Ref. [18] 

P2 = V s-^d Pj+i = 2r]Pj{r]) —Pj_i. Here we consider this parameterization for 
the gluon distribution and the corresponding Pomeron flux [Eq. (19)] as input 
in our calculations. The parameters used are from the HI fit in Ref. [18]. 



3 Results and discussion 



In the previous section, we have reviewed the formulas for the open-charm 
contribution to the proton diffractive structure functions in the perturbative 
QCD formalism and Regge based approach. In that follows, one computes 
the charm diffractive structure function considering these different analysis 
without additional parameters. In Fig. 2 one presents the results for the QCD 
approach (solid lines), the CKMT model (dashed and long-dashed lines) and 
the QCD analysis from Royon et al. [18] (dot-dashed lines) using rric = 1.5 
GeV. In particular, we consider two possibilities for the effective Pomeron 
intercept ajp{0) = 1 + A{Qljj-), which determines the xjp dependence of the 
Pomeron flux. Basically, we have considered the higher (^^^(O) = 1.24) bound 
obtained in the HERA fit and also an intercept Q^-dependent (the CKMT 
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Pomeron) . 



Regarding the CKMT approach, as the parameters have been constrained for 
— 5 GeV^, we initially compare our predictions with the experimental re- 
sults for Q'^ — A GeV^. We have that Regge based approach agrees with ZEUS 
data both in shape and overall normalization for dependent Pomeron inter- 
cept and/or fixed ajp = 1.24. For larger we have that these two choices give 
difi:erent normalizations. However, due to the scarce data a discrimination is 
still not possible. Moreover, this result can be modified by the QCD evolution, 
which is not considered in our analyzes. On the other hand, the (3 dependence 
predicted by the CKMT approach is consistent with the behavior present in 
the experimental measurements. This result is supported by the phenomeno- 
logical analyzes from ZEUS [12], which uses a fitting procedure based on QCD 
factorization for diffractive DIS in order to determine the diffractive quark and 
gluon distributions. 

The result when using the gluon distribution from Ref. [18] is quite different 
from that one coming from the CKMT model. In particular, the deviation is 
increasingly larger at small xjp and large Q^. Moreover, the behavior at small 
P has changed, which becomes fiat at this kinematical region. The reason for 
that is an almost fiat diffractive gluon distribution at small f3 coming out of 
the fit of Ref. [18], whereas one has a singular behavior when considering the 
CKMT gluon distribution on the Pomeron. These results strongly indicate 
that the charm diffractive production could allows us further constrain future 
analysis on the diffractive parton (gluon) distributions on diffractive DIS. It 
should be noticed a new set of NLO DGLAP/QCD diffractive parton distri- 
butions has been recently determined [25] (preliminary), which includes for 
the first time both experimental and model uncertainties for the error bands 
of the diffractive pdf "s. There, the small (3 behavior is steep in contrast with 
the almost fiat gluon distribution found in the fit of Ref. [18]. Therefore, it 
is expected the results using these new parton distributions will modify the 
analysis presented here. 

The perturbative QCD approach provides a steep behavior on jS in comparison 
with the Regge based one. In particular, for small (3 and small xjp the differ- 
ence between the predictions is sizeable. The main contribution in the pQCD 
approach comes from the ccg component, which is strongly dependent on the 
input for the diffractive gluon distribution. Moreover, the implicit dependence 
on f3 present in the upper limit of the /Cf-integration in Eq. (8) implies an 
additional (3 dependence. For instance, if we assume in a first approximation 
that X]pg]p{x]p,k^) oc ln/c|, it would produce a logarithmic enhancement in 
this dependence. Accordingly, the numerical calculation of Eq. (8) produces 
a strong growth at small f3 for xjp = 0.004 at both virtualities, either overes- 
timating the data points. However, the description is in agreement with data 
for xip — 0.02, even at high Q'^. It should be noticed that the QCD evolution 



11 



in the unintcgratcd ghion distribution may modify this scenario. Furthermore, 
it is important to emphasize that the pQCD approach predicts a quadratic 
dependence on xjpgjp. Consequently, for a typical powerlike behavior we ex- 
pect a stronger dependence than present in the Regge models. Therefore, a 
better discrimination between the models can be obtained by increasing the 
data statistics and enlarging the kinematical window. 

As a summary, it was shown the diffractive production of open charm is an 
important observable testing QCD dynamics. The ZEUS collaboration has re- 
cently measured [12] the open charm diffractive structure function 
which is extracted from charmed mesons D*^(2010) production. The data 
demonstrate a strong sensitivity to the diffractive parton densities. Here, we 
have contrasted the pQCD two-gluon exchange approach and Regge/QCD 
models. For the first one, the saturation model was considered in order to 
write down the unintegrated gluon distribution. In this case was observed 
good description at larger xjp, but a sizeable underestimation at smaller val- 
ues of Xjp, mostly at smaller f3, is verified. Concerning Regge approach, the 
CKMT model gives a reasonable data description in shape and normaliza- 
tion with/without a Q^-dependent Pomeron intercept. On the other hand, 
the Regge/QCD approach of Ref. [18] provides a fiat behavior at small /3, 
which is not consistent with the current measurements. A comparison with 
the most recent parameterizations of diffractive pdf's is timely. We conclude 
that an increasingly experimental statistics on this process would help to con- 
strain the diffractive gluon distribution appearing in diffractive factorization 
approaches and/or discriminate among several parameterizations for the gluon 
distribution in the Pomeron in approaches based on Regge phenomenology. 
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